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Abstract 



We explore a new definition of the persistence exponent, measuring the prob- 
ability that a spin never flips after a quench of an Ising-like model at a tem- 
perature < T < Tc, while the usual definition only makes sense at T = 0. 
This probability is now defined for spin blocks, and a general scaling for it, 
involving time and block linear size is introduced and illustrated by extensive 
simulations. 
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Recent years have seen significant progress in the study of coarsening systems [0]. Our 
understanding of phase ordering phenomena is now structured by simple ideas such as dy- 
namic scahng and universahty. For instance, it is well established that nonconserved order 
parameter dynamics are characterized by a single length scale L(t) ~ t^^^, and that a 
nontrivial exponent A appears in the scaling of the order parameter correlation function, 
{(p{'x,t')ip{'x,t)) ~ [L{t)/ L(t')]^, for t' ^ t. However, advances in this field have been to 
a large extent boosted by progress in numerical simulations of lattice systems such as the 
Ising model, while few analytical results are known except in Id. Therefore surprises are 
still to be expected as one probes more and more complicated correlations. 

Such a surprise came out recently as much interest was devoted to the study of the so- 
called persistence probability • Consider the following simple question: in a simulation 
of the Glauber dynamics of the Ising model at zero temperature, what is the fraction of spins 
p{t) which have never flipped since the initial time ? It turns out that p exhibits a nontrivial 
algebraic decay p{t) ~ t~^. A quantity such as p{t) involves the whole history of the system 
and is not easy to study analytically. Derrida et al showed analytically that 6 = 3/8 
in Id, but in higher dimensions, 6 could only be determined by numerical simulations 
or approximate methods More generally, the probability that a stochastic physical 

quantity has never changed sign since the origin of its evolution arises naturally in the 
context of nonequilibrium systems. Even for simple scalar diffusion with zero mean random 
initial conditions, a nontrivial algebraic decay is found . 

For ferromagnetic systems, until very recently, persistence had only been defined and 
studied at zero temperature for a single spin (local order parameter) p|,|^,|^J|,§] or at Tc for 
the total magnetization (global order parameter), where it yields a new independent critical 
exponent 0. For T > 0, p(t) decays exponentially due to thermal fluctuations, and the T = 
definition does not look very promising. However, there are good reasons to be interested 
in a definition at a finite temperature < T < Tc. First, numerics, renormalization group 
arguments or large- calculations assess that finite (noncritical) temperature correlations 
have the same scaling as zero temperature correlations, with the same A. Thus it is natural to 
expect the same kind of universality to hold for the persistence exponent, and it is worthwhile 
checking this point. In addition, some discrepancies were found between the value of 9 
for discrete (Ising) and continuous (y^^) models. These discrepancies were attributed to 
anisotropy effects in lattice systems at zero temperature, which should be lowered at finite 
T. Thirdly, certain important models such as the Ising model with conserved Kawasaki 
dynamics do not coarsen at zero temperature, due to finite energy barriers, and must be 
simulated at T > 0. However the corresponding continuous model (model B) exhibits 
coarsening at zero temperature, and thus one needs a rule to extrapolate information from 
finite temperature simulations using Kawasaki dynamics. 

In a very recent paper, Derrida proposed to study persistence at finite temperature 
for nonconserved Ising and Potts models by comparing two systems, A and B, evolving with 
the same thermal noise from two different initial conditions: completely random for A, and 
all spins equal to +1 for B (fundamental state). The idea is that B experiences flips solely 
due to thermal fluctuations in an ordered system. Thus a flip is recorded only when a spin 
at the same site in both samples does not flip simultaneously in A and B. Derrida found 
in d = 2 that the corresponding persistence probability p{t) decays algebraically with an 
exponent close to the value of 6* at T = 0. More extensive simulations performed by Stauffer 
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TT| also suggest a temperature-independent exponent equal to ^ in d = 2, but significantly 
different in d = 3 and d = 4. The value found in d = 3 for T > is in good agreement 
with an approximate continuous theory at T = (see conclusion). Derrida's method 
is ingenious and straightforwardly implemented, but it cannot be used to study conserved 
dynamics, as system B would not evolve with Kawasaki dynamics, and it is not easy to 
generalize to a continuous field. 

In this Letter, we propose a very natural method to study persistence at finite tempera- 
ture, namely "block scaling" , which can be directly performed on a single sample. The idea 
stems from a la Kadanoff TenoTmaAiza.tion group ideas. At finite temperature, we consider 
the persistence of coarse-grained spin variables obtained by integrating the order parameter 
(spin) on blocks. When the size / of the blocks is increased, the effective temperature flows 
to zero, which establishes a connection with the zero temperature dynamics. It is clear that 
this definition also applies to continuous models. We shall restrict ourselves to nonconserved 
ferromagnetic models (model A) with L(t) oc but the same method can be used to study 
conserved models. 

Before considering finite temperature, it is instructive to see how block scaling works 
at T = 0, for which two persistence exponents can be defined: 6 for a single spin (local 
order parameter), and 9o corresponding to the probability po(t) that the total magnetization 



(global order parameter) has never changed sign 012|. Majumdar et al [^] have shown the 



exact result = 1/4 for the Ising model in Id. Cornell and Sire performed direct 
numerical simulations of po(^) in c? = 2, by recording the time when the global magnetization 
first changes sign at T = 0. This requires a very large number of runs, which drastically 
limits the sample size (Lmax ~ 128). In addition, finite size scaling is not very conclusive, 
leading to a large uncertainty on the value of Oq ~ 0.06 ~ 0.11. 

We now show that block scaling leads to a much easier determination of 9q aX T = 0, 
before moving to finite T. Let us consider blocks of size / and the probability Pi{t) that 
the total magnetization of a block has never changed sign since t = (we will use blocks 
with odd number l'^ of spins). For large time, when L{t) 3> blocks behave as single 
spins and pi{t) ~ Qt~^, q being an increasing function of /, since obviously at large time 
Pv{t) > pi(t) if /' > /. At early times, when L(t) < the system effectively sees infinite 
blocks, and pi{t) oc t^^'\ where is the persistence exponent of the total magnetization at 
T = 0. Moreover, in the initial configuration, the larger the blocks the smaller the relative 
fluctuations of the magnetization around zero, therefore pi'{t) < pi{t), for /' > /. The cross- 
over between the two regimes should occur at t oc These remarks lead us to the following 
large / scaling, 

p^it) ~ l-'^fit/l') (0.1) 

where f{x) oc x~^° when x — and f{x) ~ when x — oo. For finite t, pi{t) must tend 
to a finite value for I oo, equal to the probability that the global magnetization never 
changed sign. This requires a = 26q. Hence, computing pi{t) for several values of / makes it 
possible to determine 6^0 by adjusting its value to obtain the best data collapse. 

To check this scaling, we simulated the T = Glauber dynamics for the Ising model 
in (i = 2 on a 2000^ lattice with blocks of linear size 1,5,9,15,19,25, and 31. 20 samples 
were averaged to obtain the final data presented in fig. |I]. We find excellent scaling, with 
6q = 0.09. Similar results were obtained in Id, confirming the scaling relation of Eq. (p.l|) 
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and the theoretical value = 1/4 (fig. 0). Therefore, block scaling is a very convenient and 
reliable method to determine 6q. 

Now let us move to a finite temperature < T < (not too close to Tc, a case studied 
in a forthcoming paper ||13|)- The difficulty in defining a persistence exponent comes from 
the fact that a spin may fiip due to thermal fiuctuations, leading to an exponential decay 
p{t) ~ exp(— t/r). Indeed, at T = 0, a spin fiips only when it is crossed by an interface 
between a + and a — domain, whereas at finite temperature, the dominating process at 
late time, when the domains are large, is the fiip of a spin within a domain due to thermal 
fiuctuations. Therefore, at low temperature, it is natural from classical kinetics intuition to 
expect an Arrhenius law r ~ exp(— A£^/T), where AS is the energy barrier to fiip a spin (or 
a block) within an ordered domain. As T — > 0, r diverges and p crosses over to a power law. 

It is instructive to justify the Arrhenius law from a random process viewpoint. Let us 
consider a block of linear size /, and spin block variables y?/. When L{t) is large enough, 
the system can be considered locally at equilibrium inside a domain, and, since there are 
no long-range correlations, the relative fiuctuation of ipi has the scaling ALpi/{ipi) oc \JtJjA. 
Thus Pi{t) is essentially the probability that a stationary random process X{t) with zero 
mean and mean square fiuctuation (X^) = T//'^ crosses a barrier of amplitude of order 1. 
If X{t) is Gaussian and Markovian, it is the solution of a simple Langevin equation, with a 
Gaussian white noise rj{t) with {rj{t)r]{t')) = 2T/l'^6{t-t'). Then it is immediatly seen that in 
exponential time u = e*, pi{u) is the survival probability of a simple Id random walker with 
diffusion coefficient 2T/l'^, starting from x = with a moving absorbing wall at x{u) oc ^/u. 
When the amplitude of the fiuctuation vanishes, i.e. for small T or large /, this survival 
probability can be evaluated by using the unperturbed solution of the diffusion equation 



At large u, piiu) decays with a power law piiu) oc and {3 oc y l'^ /T exp{—Cl''' /T) , 
where C is a constant. Thus we recover the heuristic Arrhenius law with r oc l/f3. 

The actual stochastic process (fi{t) is certainly non-Markovian. However, for / much 
bigger than the equilibrium correlation length, it is nearly Gaussian. Moreover, its correlator 
C{t) = {ipi(t)(pi{0)) can be bounded by two Markovian exponential correlators (because 
there is no long range correlation in time at equilibrium), and thus the Arrhenius law 
still holds with proper constants inserted (although the power law in the prefactor may be 
modified) [jT5[. The important point is that the effective temperature entering the Arrhenius 
law of the spin blocks is cut by a factor l'^ and that r diverges very quickly when / is 
increased, leading to a fast cross-over to the T = behavior. For t <^ r, pi{t) is expected 
to behave in the same way as for T = 0. Finally, at finite temperature (not too close to 
Tc, in the vicinity of which a different scaling arises []T3|), we expect a scaling of the form 
Pi{t) ~ 1"^^° f{t/P) exp[—t/T{l,T)], involving two cross-over times clearly visible in fig. 3, 
which shows the result of 2d simulations performed at T = 2Tc/3 on a 1000^ lattice, with 
/ = 1, 3, 5, 7, 9, 11, 13. The exponential decay is clearly visible for / = 1 and / = 3. However, 
for larger blocks, r is bigger than the simulation time, and pi{t) has the T = behavior, 
with a power law decay with exponent 6 fully compatible with the T = value {6 = 0.22), 
for t ^ and a power law decay with exponent 6*0, for t < l^, just as expected. Figure 
^ shows the scaling with 6q = 0.09 (for / = 7, 9, 11, 13, and a slightly smaller T = Tc/2 to 
eliminate the effect of the exponential cut-off). 

Thus, block scaling leads to a clear definition of 6 at finite temperature as the exponent 
of the algebraic decay of the scaling function f{x). We find that in 2d, the exponents 6 and 
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^0 do not depend on T and are equal to their T = value, in agreement with the results 
obtained with Derrida's definition. It is also very satisfactory to observe that both scaling 
functions of fig. 1 and fig. 4 are identical up to a multiplying factor, in a very similar way 
as what is known for the equal-time two-point spin correlation function 

We conclude with a look at the puzzling 3(i case. Using the present block method, we 



find an exponent 6't>o consistent with the value 9t>q = 0.26 obtained by Stauffer |]Tl| using 
Derrida's definition, but different from the T = value 9t=o = 0.17 [0,0. In fact, it is 
well-known (although a precise explanation is still lacking) |1T6|,|6[|, that the domain length 
scale L{t) does not grow as t^/^ in 3d, but as presumably due to lattice effects. If we 
now express our general scaling as a function of the more intrinsic L{t) instead of time itself, 
we find that for T = and T > (in the latter case for ^ t ^ r(/, T)), both persistence 
probabilities decay as p{t) ~ L(t)-^ with the same 6 ^ 0.17/0.33 ^ 0.26/0.5 ^ 0.52, in 
good agreement with the theoretical prediction of @]. 

We are very grateful to B. Derrida and S. Cornell for helpful discussions. 
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FIGURES 



-0.6 
"f -0.7 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

log{t/l') 

FIG. 1. Block persistence at T = obtained from simulation of the nonconserved Ising model 
on a 2000^ lattice, for Z = 1, 5, 9, 15, 19, 25, and 31 (from bottom to top in the insert). pi{t) decays 
as t~^° at early time and as at large time. Excellent scaling is then obtained taking = 0.09. 
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FIG. 2. Similar scaling as in fig. 1, for a Id spin chain (200000 spins, 10 samples), with block 
size Z = 1, 21, 41, 61, 91 (from bottom to top in the right part of the insert). Z = 1 is omitted in the 
scaling, and the data collapse improves as the block size increases. 
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FIG. 3. pi{t) for T = 2Tc/3, and block sizes / = 1, 3, 5, 7, 9, 11, 13. 
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FIG. 4. pi{t) expressed in scaling form for T = Tc/2, and block sizes / = 7, 9, 11, 13, using the 
same value for 29o = 0.18 as in the T = case. Note the similarity with the T = scaling function 
of fig. 1. 
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